Nonlinear ac conductivity of interacting Id electron systems 
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We consider low energy charge transport in one-dimensional (Id) electron systems with short range 
interactions under the influence of a random potential. Combining RG and instanton methods, we 
calculate the nonlinear ac conductivity and discuss the crossover between the nonanalytic field 
dependence of the electric current at zero frequency and the linear ac conductivity at small electric 
fields and finite frequency. 
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I. INTRODUCTION 

In Id electron systems, the effect of both interactions 
and random potentials is very pronounced, and a vari- 
ety of unusual phenomena can be observed^'^. The lin- 
ear dc conductivity shows a power law dependence on 
temperature T at higher temperatures'^''* , but is ex- 
ponentially small at low temperatures and vanishes at 
zero temperature^'®. The ac conductivity vanishes like^ 
~ a;^(ln(l/w))^ and shows several cross-overs to other 
power laws at higher frequencies^. 

Much less is known about the non-linear conductivity. 
At zero temperature and frequency charge transport is 
only possible by tunneling of charge carriers, which can 
be described by instanton formation. The nonlinear dc- 
conductivity is characterized by^"^^ / ~ exp(— y^o/^E) 
provided the system is coupled to a dissipative bath. 
Without such a coupling, the current was recently sug- 
gested to vanish below a critical temperature"'^'^'^'*. 

In this work, we calculate the low energy non- linear ac 
conductivity for systems with random pinning potentials 
and discuss the crossover between linear ac response at 
small fields and nonlinear dc response at large fields. To 
be specific, we consider a charge density wave (CDW) 
or spinless Luttinger liquid (LL) pinned by a random 
lattice potential which can be described by the quantum 
sine-Gordon model with random phases. We first scale 
the system to its correlation length, where the influence 
of the potential is strong and a semiclassical instanton 
calculation becomes possible. The response is dominated 
by energetically low lying two level systems (TLS), whose 
dynamics is described by a Bloch equation. 

Microfabrication of quantum wires or Id CDW 
systems^^ should allow to test our predictions experimen- 
tally. Indeed there is a number of recent experiments on 
carbon nanotubes^®^*^ and polydiacetylen^^ which seem 
to confirm the variable range hopping prediction for the 
dc-conductivity made in^*''^*. 



II. AC CONDUCTIVITY OF ID DISORDERED 
SYSTEMS 

In the following, we present a heuristic derivation^ of 
the Mott-Halperin result^ for the ac conductivity of a 



one-dimensional disordered electron system without in- 
teractions. In the end of the section, we indicate how 
this result can be generalized to interacting electrons. 

In one spatial dimension, all electron states are local- 
ized and wave function envelopes decay on the scale of 
the localization length ^loc. We divide the system into 
segments with size ^/oc. The typical energy separation 
of states within one segment is the mean level spacing 
A — 1/(pf6oc), where is the density of states at the 
Fermi level per unit length. Levels in neighboring seg- 
ments are coupled by the Thouless energy t(^ioc) = A. 
When we consider the coupling between more distant 
segments of separation L the coupling is reduced to 
t{L) = A exp(— L/^/oc). The coupling splits (almost) de- 
generate energy levels in different segments by an amount 
AE = 2t[L). Due to the couphng, the eigenstates of 
the Hamiltonian are even and odd linear combinations of 
states localized in segments a distance L apart. A spa- 
tially constant ac electric field E{t) = EqCOSujI causes 
transitions between levels with a separation AE = Huj, 
hence we demand 2t{L) = huj and therefore 



Lxiuj) = C;oc ln(2 A/ft w) 



(2.1) 



According to Fermi's golden rule, the transition rate for 
exciting the even linear combination to the odd one is 



given by l/r = ^ \elE'^\x^^,+nuj\'^ Pf S.ioc- Here, the 
matrix element Xt^tJ[-hu> of the position operator has to 
be calculated between the ground state with energy e 
and the excited state with energy e + Huj, it is found to 
be equal to the spatial separation Lx{oj) of the two local- 
ized states. When calculating the rate of energy absorp- 
tion i(Tac(w)£^o due to the excitations of such two level 
systems, one has to take into account that each photon 
carries the energy ftw, that only transitions from unoc- 
cupied states to occupied states are possible, and that 
the occupation probability for a state with energy e is 
determined by the Fermi function /(e). In this way, one 
finds 

/2tt 



i- ^loc- 
ft 



(2.2) 



This result can be generalized to an interacting elec- 
tron system^® by remembering the basic idea of bosoniza- 
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tion: the charge density is defined as the derivative of a 
displacement field, and a localized electronic state corre- 
sponds to a localized kink in the displacement field. In 
addition, the density of states pp at the Fermi level has 
to be replaced by the compressibility « = |^ . With these 
modifications, the above derivation can be repeated and 
one obtains an result analogous to Eq. (2.2). 

It is worthwhile to remark that the ac-conductivity 
(2.2) can be rewritten as 



(2.3) 



where ^ is the correlation length of the system. This 
result resembles the form of the result for the non-linear 
rfc-conductivity 

c^..{E) ~ .0 ^ e-^^(^)/«- (b^^" . (2.4) 



Here Lx{E) = ^loc 



denotes the spatial dis- 



eoSoiioc 

tance of the energy levels between which the tunnel- 
ing events take place, it follows from a variational 
treatment^' The prefactor fiwe-ph/eoi^o^oc takes 
into account that the dissipation rate is controlled 
by the typical frequency for electron phonon coupling 
We-ph and that the current is hence of the order^^ ~ 
eoWe-phC"^^"^^^^^'"" . We note that depending on the 
details of the electron phonon coupling, We-ph may de- 
pend on the external electric field. 

Finally, at finite temperatures, Mott variable range 
hopping gives a linear lic-condTictivity which follows from 
(2.4) by replacing L^{E) by L,(T) = ^A/ksT, i.e.''^° 



(Tdc{T) 



(2.5) 



The cross-over between the three expressions (2.3)-(2.5) 
can most easily be understood by the dominance of a 
shortest tunneling distance, as will be discussed further 
below. 



III. THE MODEL 

The models we analyze are defined by the euclidean 
action 



dxj \dy 



(3.1) 



-2u cos{p(p + 2tt(^{x)) + '^^^^ ipE{y) 



nv 



where we have rcscalcd time according to vt y, and 
(3 = l/ksT. The dissipative part of the action describes 
a weak coupling of the electron system to a dissipative 
bath, for example phonons. It is needed for energy re- 
laxation in variable range hopping processes and for 



equilibration in the presence of a strong ac field. We as- 
sume it to be so small that it does not influence the RG 
equations for the other model parameters significantly. 
The smooth part of the density is given by ■^dx'f, and 
p = 1,2 for CDWs and LLs, respectively. We consider a 
CDW or LL with Gaussian disorder, which is described 
by C,{x) equally distributed in the interval [0, 1] with cor- 
relation length equal to the lattice spacing a. 

For K > Kc{u) the potential is RG irrelevant and de- 
cays under the RG flow, while for K < Kc{u) the po- 
tential is relevant and grows, here^ Kc{Q) = 6/p^ ■ We 
assume K < Kc{u) and scale the system to a length 
^ = ae' , on which the potential is strong. After the 
scaling process, the parameters K, v, and u in Eq. (3.1) 
are replaced by the effective, i.e. renormalized but not 
rescaled, parameters K^g, t^eff, and Ues- 

We note that the ratio K/v and hence the compress- 
ibility K = is not renormalized due to a statistical 
tilt symmetry^ ^. The compressibility k = |^ is used as a 
generalized density of states for interacting systems. Our 
calculations arc valid for energies below the generalized 
mean level spacing Aq = 

In this RG calculation, we do not attempt to treat 
a possible nonlinear dependence of coupling parameters 
on the external electric field. The full inclusion of the 
external field in an equilibrium theory is not possible 
as it renders the groimd state of the system unstable. 
The quantum sine-Gordon model has an infinite number 
of ground states connected by a shift of the phase field 
by 2tt/p. Here, we concentrate on renormalizing each 
of these ground states separately and take into account 
the coupling between different ground states due to the 
external electric field in the framework of an instanton 
approach. 



IV. DISORDERED LL OR CDW 



The wall width 1 / ^/jPu^ w ^ of an instanton solution 
to the action Eq. (3.1) is for weak external fields much 
smaller than the extension of the instanton. Hence, the 
instanton action can be expressed in terms of the do- 
main wall position X{y). The discussion of instantons 
in the case of random pinning is more involved than e.g. 
for periodic pinning^^ and the calculation of closed form 
instanton solutions is not possible. For this reason, we 
look for approximate instanton solutions with a rectan- 
gular shape and extensions L^, Ly in x- and y-direction, 
respectively. As the disorder is correlated in time but not 
in space, instanton walls in x and y direction contribute 

T — ) LxSy and LySx{x) to the action, respectively. While the 
n c . . 2-K ■ " , , 



surface tension Sy = 



surface tension 



-mp — is essentially constant, the 

has a strong and random position de- 



pendence. To calculate the statistical properties of s^, we 
make use of the exact solution^^ of the classical ground 
state of a LL or CDW with random pinning in the fol- 
lowing. 
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In the limit /Toff ^ 1, quantum fluctuations arc 
strongly suppressed and the (classical) ground state of 
the model Eq. (3.1) can be determined exactly^". After 
renormalization to the scale the effective action can be 
rewritten as a discrete model on a lattice with grid size 
^, and the integration over x can be replaced by a sum- 
mation over discrete lattice sites i = x/^. In the classical 
ground state, the solution (p{x, y) does not depend on y 
any more and the y integral in Eq. (3.1) simply yields 
an overall factor v^sW- Dividing the action by /i/3, one 
obtains the classical Hamiltonian^^ 



Lx X Ly-instanton^° 



Ao 
27r2 



X] [{^i+l-^if - 2^^Weff C0s(p<^i-27rCi)] • 
i=l 

(4.1) 

Here, Woff is the disorder strength with UcS ^ 1 and 
Ci G [0, 1] is a random phase. In the effective Hamilto- 
nian Eq. (4.1), the disorder term dominates the kinetic 
term and the classical ground state of the system can be 
explicitly constructed^^. One minimizes the cosine po- 
tential for each lattice site by letting p(pi = 2Tr(Q + n^) 
with integer n?. The set of integers is chosen in such 
a way that the elastic term in Eq. (4.1) is minimized, 



AS 

IT 



= (Sx(io) + Sx(«o + L^)) Ly + 2syLx, s^(i) = Sygi. 

(4.4) 

In a noninteracting electron system, a pair of sites with 
5i <C 1 corresponds to an electron with energy — ei = 
—giAo just below the Fermi level at position io, which 
can hop to an unoccupied level with energy 62 = 52 Aq 
just above the Fermi level at position io + Lx- The trans- 
lation between the language of noninteracting elctrons 
and the bosonic language used in this calculation is sum- 
marized in Table I. 

Typically, the two lowest gi in an interval of length Lx 
are of the order l/L^, and the boundaries of a typical 
instantons will be at positions with a small surface ten- 
sion Sx ~ SyS^/Lx- Taking into accoimt the contribution 
of a dc external electric field, the total action of a typical 
instanton is 



SiLxi Ly) 



= 2Sy— Ly -\- ISyLx —L^Ly . (4.5) 



1 



Extremizing the action with respect to L^, Ly one finds^° 



: m 



i<j 



+1 



(4.2) 



Here, [(]g denotes the closest integer to (, and m is 
an integer parameterizing the infinitely many equivalent 
ground states. Excitations of the ground state change 
n° ^ n° ± 1 for sites with io < i <io + L^/^, they bifur- 
cate from one ground state characterized by m = mo to 
another ground state with m = ?no±l. The potential en- 
ergy necessary for a bifurcation at position i is according 
to Eq. (4.1) 



AH(i) = ||(2Vp)'{l±2((Ci 

= hVeSSygi 



Ci-i)-[Ci-o-iy} 

(4.3) 



with a random gi G [0,2]. Defining the localization length 
6oc = 2^|f^, one has Sy = ^. 

Quantum effects are due to the time derivative in the 
action Eq. (3.1) and give rise to tunneling between the 
ground state and excited states in the presence of an 
external electric field. Similar to the case of periodic 
pinning^^, these tunneling events are described by in- 
stantons. In the following, we describe how the action of 
a quadratic instanton can be calculated. 

The action of a bifurcation with extension Ly is just 
AH{io)Ly/vee- The action of a wall with constant y 
and length Lx can be calculated by an analogous con- 
sideration if one introduces a lattice of grid size ^ in y- 
direction. As the disorder is correlated in time direction, 
one needs not consider random phases Ci and finds an 
action hsyLx- Adding up the contributions from all four 
walls of an instanton, one finds the action of a rectangular 



LxiE) = 



2tt 



Ln — 



pneoEo pk^bqEq 
The creation rate of these instantons is 



(4.6) 



random 



-2L,(£;o)/6o 



(4.7) 



leading to the result Eq. (2.4) for the conductivity. 



TABLE I: Translation between boson and electron language. 



electron language 



boson language 



particle 

hole 
p{x) 

surface tension ^^.('to) 
surface teusicni .s,r(/() + L,, ) 



kink 
antikink 

energy -ei(io) 
eiK-rgy f_.(/ij + Ly 



Next we consider an ac field E{t), which upon ana- 
lytical continuation it ^ t turns into a field E{t). In 
imaginary time, the electric field has to obey the same 
periodic boundary condition E{t + P) = E{t) as other 
bosonic fields, e.g. the displacement field ^pij). This 
boundary condition is respected by a discrete Fourier 
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representation 



£(t) =rVi;(u;„)e 



with Matsubara frequencies a;„. A monochromatic ex- 
ternal field is hence described by E{y) = Eq cosojny, 
where time is rescaled as y = WeffT and frequency as 
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FIG. 1: Hopping of kinks from position xo to xq + Lx or back 
(full lines) in a time interval of length Ly. The lines parallel to 
the i/-axis do not contribute to the action if the line tensions 
Sx are neglected. 



LJn = tOn/vcS- In the end of our calculation, we ana- 
lytically continue Matsubara frequencies to retarded real 
frequencies iuJnVcS lu + irj. 



Which type of instantons determines the current in 
the presence of an ac field with period = 27r/tZ;„? 
In the limit of a very small external electric field, the 
only external length scale in the problem is the ac pe- 
riod L^, and the contribution of typical instantons to 
the linear response can be estimated by assuming that 
Ly has to be of order L^^. A typical instanton ob- 
tained from minimizing Eq. (4.5) with respect to for 



a fixed Ly = L^^ and for vanishing Eq has L^ 
For this solution, the current would be proportional 
to exp[— |f (1 -I- Ao/huj] and vanish nonanalytically 
for small frequencies. However, according to the Mott- 
Halperin law, the true frequency dependence should be 
proportional to^ uj'^ In^ {I / uj) . We conclude that typi- 
cal instantons do not yield the leading contribution to 
the current and that a discussion of rare instantons is 
needed. 

Indeed, besides typical instantons with w Sy£^/Lx, 
there are rare instantons with an exceptionally low 
Sx{i) + Sx{i + Lx). Such a pair of sites i and i + L^ allows 
for the hopping of a kink without changing the kink's 
potential energy much. The potential energy difference 
ei + £2 between two sites can become arbitrarily small 
in sufficiently large samples. For the following consider- 
ations, we will set it to zero in the sense that it is much 
smaller than any other energy scale in the system. In the 
discussion of the dc electric field, quantum fluctuations, 
i.e. spontaneous creation of typical instantons in the ab- 
sence of an external field, were unimportant. For pairs 
of sites with exceptionally low surface tensions, quantum 
fluctuations are important and have to be taken into ac- 
count. Here, the spontaneous formation of instantons 
describes the physics of level repulsion^'*. In our approx- 



imation of vanishing Sx, the instanton action does not 
depend on the extension Ly in time direction any more, 
hence the occurrence of single domain walls of length Lx 
with constant y is possible. Such a domain wall describes 
the hopping of a kink across the distance Lx, and its ac- 
tion is 



S single I ^ SyLx 



(4.9) 



To obtain the partition function for this tunneling de- 
gree of freedom, we must sum over all possible domain 
wall configurations in the interval L^^. A configuration 
with three hopping events is displayed in Fig. 1. Sum- 
mation over all possible configurations yields 



ZiL^) 




^-nSyL^ (4.10) 



The integration measures in the y^-integrals are normal- 
ized by ^-R'off rather than by ^ because the short time 
cutoff ft/Ao is determined by the high energy cutoff Aq 
and not by the short distance cutoff ^. We note that the 
exponent of the outer exponential function is positive, 
indicating a lowering of the ground state energy due to 
frequent tunneling between the degenerate states. The 
summation over all possible instanton configurations de- 
scribes the quantum mechanical effect of level repulsion. 
Coupled energy levels repel each other and are separated 
at least by 

h{Lx) = f^\nZ{L^) = Aoexp(-L,/eoc) ■ (4.11) 

The probability to have exactly one instanton (hopping 
of a kink forth and back) within the time interval L^^ is 
given by 



p,{Lx) = e-^^«^==/Z(L^) 
L 



(4.12) 



exp 



e ^ " 



'ZSyLx 



The optimal length Lx for such an instanton is found by 
minimizing the exponent in Eq. (4.12) with respect to 
the tunneling length. We find 



Lx — — hi ■ 



off 



(4.13) 



Using this expression for Lx, we find the probability for 
having exactly one instanton of length 



_2 /2e?^eff 



(4.14) 



5 



This proportionality of the tunncUng probabihty to fre- 
quency squared is the essence of the Mott Halperin con- 
ductivity Eq. (2.2). 

With the knowledge of the probability Eq. (4.14) for an 
instanton in resonance with the external field, we can now 
set up a calculation of the ac current. It is calculated as 
a derivative I{x,ujn) = ^^ sa{x^-u ) partition 
function with respect to the vector potential a{x, w„) = 
E{x,uin)/'^n- The field (f) couples to the vector potential 
via 

SE/h= ^ [ dx -J— ^a{x, -U)n){-U)n)(t>{x,U)n) , 
TTfi J pVeS ^ 

(4.15) 

where /3 denotes the inverse temperature. Hence, the 
current is given by 

I{x,Cb^) = -eo QnVeS jDi^^,-sw • (4-16) 

In the low energy regime, we do not perform the full 
functional integral over in order to evaluate the par- 
tition function. Instead, we sum over the relevant tun- 
neling degrees of freedom. We label such a degree of 
freedom by the position «o of the first weak link and by 
the distance between the first and the second weak 
link. The field <j){io) = + 1) = ••• = '/>(«o + L^) 
takes the values 00 and 00 + 27r/p, respectively. Further- 
more, we make use of the self averaging properties of the 
current. Instead of averaging it over the position x in 
the system and let all different types of TLSs contribute 
to it, we calculate the contribution of one TLS and av- 
erage over the parameters L^, Sx{io), and Sx{io + L^)- 
So far we considered instantons with vanishing surface 
tension Sx = 0. We now extend these considerations 
to instantons with surface energies smaller than hio, i.e. 
we are concentrating on sites with gi < Keff£,(jJn, where 

(«o) = (io + ix) = Syg2. The probability that 

the position .x, for which we want to evaluate the current, 
is inside an active instanton is L^;/^ times the probability 
for finding a weak tunneling link at a given site. In this 
way, we obtain for the average current 



(^X'^n) 



dyi dy2 



M<^n;yi,y2) e-^«[^°l(4.i7) 



Here, we evaluate the current under the approximation 
that there is exactly one instanton (two tunneling events 
at times yi and 7/2) in the interval of length L^. The two 
gi integrals contribute a factor (Kee^i^n)^ as we neglect 
the dependence of Se on the g^. The integral over is 
evaluated by the saddle point method just taking into ac- 
count instantons of the optimal length L^; = -!- In . 

Sy jit^J\eti 

As pi is a function of Lx/^ioc, we have to use the same 

variable in the integration measure to perform a saddle 
point approximation and obtain a factor ^loc/C from this 



transformation. If yi is the time of the first tunneling 
event and y2 the time of the second tunneling event, we 
define the new variables y = (yi+?/2)/2 and Ly = 2/2 — J/i- 
The Fourier transform of the displacement field (j) for such 
an instanton is given by 



cj>{u,n;Ly,y) = -e'---y sin^ 



(4.18) 



For this field configuration, the coupling to the external 
electric field contributes the action 

SE{Ly,y)/h = sm{u)ny) sin(d;„Ly/2) . 



(4.19) 



The barrier size is not fixed by the electric field as in 
the dc limit, and one has to consider both forward and 
backward jumps as for the standard thermally assisted 
flux flow argument^^. Hence, the probability for the in- 
stanton to be in phase with the external field is given by 
2smh[SE{Lu;, Lx)/h]. Expanding to linear order in Se 
and integrating over Ly and y, we find for the current 



-eoVeSi^n 



U)n^KeS^2 



dy 



8 e, 



"dLy 2SEiLy,y) 
^(u}n;Ly,y) 

(4.20) 



After the analytical continuation a)„Ueff — * —ii^ the real 
part of the conductivity agrees with formula Eq. (2.2) 
up to a numerical factor. As both ^loc and are 
proportional to Kgg, the conductivity contains a factor 
of K^g in agreement with the result 



8,26 



V. DESCRIPTION BY BLOCK EQUATION 

The instanton calculation presented in the last section 
needs to be improved upon in two respects: first, the 
classical level separation parameterized by the gi was not 
fully taken into account, and second, nonlinear correction 
in the strength of the external field are not considered yet. 
To achieve these goals, we use the concepts developed in 
the last section in a real time quantum mechanical cal- 
culation. Instantons in the imaginary time formalism 
correspond to the hopping of kinks from one level just 
below the chemical potential to another level just above 
the chemical potential. The localized states of a Id dis- 
ordered system can be modeled by an ensemble of these 
TLSs, and the average properties of the system can be 
calculated by averaging over the parameters of the TLSs. 

We consider a TLS with spatial extension and on 
site energies — ei and £2 with = 5,Ao. In a system 
of noninteracting electrons, the negative energy — ei cor- 
responds to a particle below the Fermi level, and the 
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positive energy £2 to an unoccupied site above the Fermi 
level. The two sites are coupled by a distance dependent 
hopping integral to(-^x) according to Eq. (4.11). Such a 
TLS is described by the Hamiltonian H = Hq + He with 

-^0 = ^{ii+e2)(Tz-to{Lx) (Jx,He = ^ eoLxEo COS utaz ■ 

(5.1) 

The position of the tunneling kink is measured by the 
operator x = ^L^cTz, and the current operator is given 
by -f = p(ei, £2, Lx)eox. Here, 



1 dei de2 dL^ 



^AoAo e 



(5.2) 



denotes the spatial density of TLSs with given parameter 
values. Ho is diagonalized by the unitary transformation 



H — > exp{iipay/2)H exp{—i(pay/2) 



(5.3) 



with ip = arctan In the new basis, Hq corresponds 

to a static field in ^;-direction, and He to an oscillat- 
ing field with x component proportional to sin (p and z- 
component proportional to cos tp. 

In principle, the transformed Hamiltonian should now 
be solved in a noncquilibrium setup in a dissipativc envi- 
ronment. In general, this type of problem is difficult to 
deal with in full generality^''. However, the problem sim- 
plifies if one does not treat an individual quantum system 
but averages over a whole ensemble instead. Such an en- 
semble of TLSs or spins interacting with an oscillatory 
electric field and subject to relaxation processes can be 
described by Bloch equations^^. We denote the ensemble 
polarization of the TLSs in the transformed basis by the 
pseudospin vector p. The current / is then proportional 
to the the pseudospin component in ^/-direction, 

I = -^PyPi^i, £2, Lx)eoLj-^^^j^ . (5.4) 

The pseudospin polarization p of the TLSs follows the 
Bloch equation 

^^p = -— {p-p^ + apx E . (5.5) 
at — tq — — 



Here, E = 



Eq sm{ip) cos{ijjt), 0, -|- Eq cos{ip) cos{u)t) 
V(e^ + ei)/4 + t^(ix), and 



[0,0,-1], A 



^0 

a = eoLx/h. Inelastic processes are described by the 
phenomenological damping constant tq. In the absence 
of an external electric field, the pseudospin relaxes due 
to this damping to its equilibrium value p^, i.e. the 
particle is in a superposition of states localized at to and 
ia + Lx- At temperatures much lower than the hopping 
integral to{Lx), the dissipative bath cannot destroy this 
coherent superposition and localize the particle, as the 
localized states have a higher energy than the symmetric 
linear combination. 



The solution of Eq. (5.5) is described in detail in 
reference^* and we do not reproduce it here. We find 
that to order 0{Eq), one type of TLS contributes to the 
conductivity 



CTTLs(w;ei,e2,ix) = P^i,^2,L 
1 - 



w — 



2A 



elLlt\Lx) 
1 {eoLxEoto{Lx)f 



?i2A2 (^_^)2+ 1 



.(5.6) 



The reduction of the linear conductivity becomes effec- 
tive for strong ac fields, when both states of the TLS 
are occupied with comparable probability. In order to 
calculate the conductivity of the disordered sample, we 
integrate over all possible; parameter values ei, £2, and 
Lx and obtain the final result 

ReaM = crolLl{uj)i!u.nni-2^^^0^) (5.7) 

with the optimal tunneling length given by 

2Ao 



i'x(w) = 60c In 



(5.8) 



The linear part of Eq. (5.7) is proportional to K^g and 
agrees with the result of Fogler^ . This linear conductiv- 
ity describes the response of a disordered Id system in 
region MH of Fig. 2. For an unscreened Coulomb inter- 
action, in Eq. (5.7) one factor of hut has to be replaced 
by^^ eQ/eLx{co), where e is the dielectric constant of the 
system. 

When eoEoLx{uj) « h/ro, higher order terms become 
important and the ac current will saturate as a function 
of £^0- The value Eg of the electric field where the current 
saturates can be estimated from Eq. (5.7) as 



E.. = 



h 



eQToLx{ijj) 



(5.9) 



As the nonlinear conductivity is defined as the ratio of 
current and electric field, in the saturation regime one 
obtains 



P 



C7ac(w, E> Es)= ao-LxiE) 



S, tqAo 

(5.10) 

The region in w-E space where the nonlinear conductivity 
Eq. (5.10) can be observed is labeled MHS in Fig. 2. 



VI. DISCUSSION 

How does the linear conductivity Eq. (5.7) connect to 
the creep current in strong fields? The calculation of 

the nonlinear dc conductivity involves the optimal length 
scale Lx{E) in Eq. (4.6) for tunneling processes. The 
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FIG. 2; Different transport regimes in a disordered CDW or 
LL with dissipation as a function of field strength in units of 
£^throshoid ~ Afl/peo A'off ^loc and frequency: adiabatic quan- 
tum creep in A, quantum creep with increased current noise 
in QC, ac conductivity following a Mott-Halperin law in MH, 
Mott-Halperin law with saturation for large field strengths in 
MHS. 



crossover from ac to dc conductivity takes place when 
the two length scales L^iuj) and Lx{E) match, i.e. for a 
crossover frequency 



^„-ix(Bo)/«io 



(6.1) 



For id ~ iOx ; the magnitude of the creep current /creep ^ 
exp(— 2La;(i?)/^ioc) agrees with the magnitude of the ac 
current described by the conductivity Eq. (5.7), as the uj"^ 
term in Eq. (5.7) matches the exponential dependence on 
field strength of /creep- Then, the expression Eq. (5.10) 
turns into 



this estimate agrees with a more sophisticated calcula- 
tion, in which a TLS is coupled to phonons with an 
Ohmic spectral function'^". 

In the crossover region lo ~ lox between regimes QC 
and MHS in Fig. 2, there are two different types of 
TLS contributing to the current. In the ground state 
of a typical TLS with bare energy difference ei + 62 ~ 
^0^/ Lx{E) ^ t{Lx), most of the charge is localized at 
one of the levels. Under the influence of an external field, 
the charge hops irreversibly from one level to the other , 
as the hop is generally accompanied by an inelastic pro- 
cess. On the other hand, the ground state of a TLS with 
exceptionally low bare energy separation ei + £2 5: i{Lx) 
is the even parity combination of wave functions centered 
around the individual levels, and absorption of a photon 
excites the TLS to the odd parity state. 

In the quantum creep regime, the time dependence of 
the current is calculated using the time dependent field 
in the formula for the dc current /creep- The adiabatic 
regime (region A in Fig. 2) is reached for frequencies 
smaller than the dc hopping rate /crccp/eo- While the 
average current is stronger than the current noise in the 
adiabatic regime, the current noise is stronger than the 
current for /crccp/eo < w < Wx. 

In summary, we have discussed the crossover from a 
nonlinear creep current in a static electric field to the 
linear ac response in Id disordered interacting electron 
systems. While the linear ac conductivity is described 
by a generalized Mott-Halperin law, for stronger fields 
one finds a reduction of this linear conductivity as both 
states of a TLS are occupied with comparable probability. 
The crossover between nonlinear ac conductivity and dc 
creep current occurs when the spatial extension of TLSs 
matches the length scale for tunneling of kinks. 



8roAo e 



-2L,(£;)/Cio 



(6.2) 



providing us with an estimate of the prefactor of the ex- 
ponential factor describing dc creep. Identifying H/tqAq 
as a dimensionless measure for the dissipation strength. 
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